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Abstract. The charge-carrier transport properties of ultrathin metallic tilms are analysed with ab-initio methods using 
the density functional theory (DFT) on free-standing single crystalline slabs in the thickness range between 1 and 8 
monolayers and compared with experiments for Pb films on Si(lll). A strong interplay between bandstructure, quan- 
tised in the direction normal to the ultrathin film, charge-carrier scattering mechanisms and magnetoconduction was 
found. Based on the bandstructure obtained from the DFT, we used standard Boltzmann transport theory in two dimen- 
sions to obtain results for the electronic transport properties of 2 to 8 monolayers thick Pb(l 1 1) slabs with and without 
magnetic field. Comparison of calculations and experiment for the thickness dependence of the dc conductivity shows 
that the dominant scattering mechanism of electrons is diffuse elastic interface scattering for which the assumption 
of identical scattering times for all subbands and directions, used in this paper, is a good approximation. Within this 
model we can explain the thickness dependences of the electric conductivity and of the Hall coefficient as well as the 
anomalous behaviour of the first Pb layer. 

PACS. 73.50.Jt Electronic transport phenomena in thin films: Galvanomagnetic and other magnetotransport effects - 
73.61.At Electrical properties of specific thin films: Metals and metallic alloys - 73.20.At Electron states at surfaces 
and interfaces -71.15 .Mb Density functional theory 



1 Introduction 

Spatial confinement of the electron motion in the direction nor- 
mal to the film causes quantisation of the perpendicular mo- 
mentum and consequently a discrete set of electron energy sub- 
bands. This leads to a variety of quantum-size effects (QSE), 
like oscillation of the charge-carrier density, Fermi energy, elec- 
trical conductivity and of the Hall coefficient as a function of 
film thickness, to mention only a few of them. Although ap- 
preciable progress has been made in understanding the QSE, 
many aspects are still unclear. Here we will discuss the elec- 
tronic band structure and its consequences on the electronic 
charge transport with and without magnetic field in ultrathin 
films. 

The classical size effect on the conductivity of thin films is 
explained with the Fuchs and Sondheimer theories [1-3] which 
relate the electrical conductivity in the film and surface or in- 
terface scattering. Their approach is based on the Boltzmann 
transport equation for bulk electrons that are partially reflected 
and partially scattered on the film surfaces. Assuming diffuse 
and elastic surface scattering they find a linear thickness depen- 
dence of the conductivity, often in agreement with experiments 
even for film thicknesses d much smaller than the electron elas- 
tic mean free path £, although the picture of 3-dimensional elec- 
tron motion with specular reflections at the surface or interface 
makes sense only when d >> £. A linear dependence of the 
conductivity has also been obtained by Camblong and Levy 



[4], who used the Kubo formalism and assumed a free-electron 
model with diffuse scattering on the boundaries. 

A quantum size effect is the oscillation of the conductivity 
with thickness d, observed in the layer-by-layer growth mode 
of thin films. These oscillations have the period of the per- 
pendicular lattice spacing and are caused by the roughness of 
films with partial monolayers that are not closed. For film thick- 
nesses d « I the discrete nature of the electron energy lev- 
els becomes essential and other quantum-size induced effects 
emerge on top of the layer-by-layer induced oscillations. The 
transport of free electrons in confined geometries has been dis- 
cussed by Sandomirskii [5] and later by Trivedi and Ashcroft 
[6]. They predicted a saw-tooth like oscillations of the in-plane 
conductivity with thickness. The period of the oscillation is 
equal to half the Fermi wavelength Xp [6] . In reaUty, however, 
the film thickness increases in monolayer steps and not contin- 
uously as assumed in the above models. Therefore, one expects 
a range of thicknesses with large QSE amplitudes when d is 
close to a multiple of Ai?/2, separated by regions with small 
QSE. 

Apart from the theoretical attempts mentioned above, there 
have been several others which describe charge transport in thin 
films [7-13]. Common to most of them is the assumption of 
free electrons. However, models based on free electrons break 
down if they are used to explain the galvanomagnetic phenom- 
ena (like the Hall effect), which depend on details of the band 
structure and on the electron/hole ratio, as will be shown below. 
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More quantitative models are needed that take into account the 
detailed band structure of the film under size quantisation. 

In this paper we address the general issue of charge-carrier 
scattering, transport and Hall effect in ultrathin metallic films. 
The Hall coefficient has the advantage that it is more sensi- 
tive to the QSE-mediated changes in the bandstructure and less 
sensitive to the details of the scattering mechanisms than the 
electrical or Hall conductivities. The model calculations will be 
performed on free-standing Pb(l 1 1) slabs and compjired with 
experiments for Pb on Si(lll). As substrates, semiconductors 
are used because they have no electron states at the Fermi en- 
ergy to which the metal conduction electrons, responsible for 
the charge transport, could couple. Semiconductors are also 
good insulators in the relevant temperature range (up to 150 K) 
so that they can be considered as an insulating substrate. De- 
pending on substrate and annealing conditions, both crystalline 
and semi-amorphous films can be grown. This is important be- 
cause it allows a comparison of magnetoconductive properties 
of crystalline and non-crystalline films. For the Hall coefficient 
the degree of crystallinity is of minor importance. 

The bandstructure of free-standing ultrathin Pb(ll 1) films 
has been calculated by Saalfrank [14] and more recently by 
Matrezanini et al. [15]. They discussed the role, QSE has on 
the oscillations of the Fermi energy ep and of the density of 
electron states np at cp- but did not consider the charge-carrier 
transport phenomena. Charge transport properties at low tem- 
peratures are dominated by scattering at bulk and surface im- 
perfections. For most ultrathin metallic films with d < £, scat- 
tering at the surface roughness is much stronger than scatter- 
ing at bulk imperfections, the electrical conductivity is limited 
by roughness of the boundaries to the vacuum and to the sub- 
strate. The variation of the conductivity with thickness d of ul- 
trathin rough metallic films has been investigated by Fishman 
and Calecki [9, 10]. In case of free electrons and uncorrelated 
surface height fluctuations, they obtain a rf^'^ thickness depen- 
dence of the "residual" conductivity a = [p{d) ~ p(oo)]~^, 
where p{d) is the in-plane film resistivity, and a Hall coeffi- 
cient Rh proportional to d. None of these dependences agrees 
with the experiments on ultrathin Pb films on Si(l 1 1), where an 
approximately linear thickness dependence of a and a strongly 
oscillating thickness dependence of Rh has been reported [16- 
18]. 



2 The Model 

We treat metallic films as free-standing, ideally ordered slabs. 
The density functional theory (DFT) is applied to calculate the 
electronic band structure, density of states and Fermi lines. Us- 
ing the Boltzmann transport equation in the relaxation time ap- 
proximation and in two dimensions we then calculate the trans- 
port properties of ultrathin slabs. We will see that most of the 
charge-carrier transport properties as a function of layer thick- 
ness d are already reproduced in such approach. 

The basis of our charge-transport calculations are the elec- 
tron band structures. In a slab, the translational symmetry in the 
direction normal to the slab is broken, k± is no longer a good 
quantum number, and each band splits into d discrete subbands. 
One can imagine these subband states as standing waves be- 
tween the two surfaces. Of course, the subbands transform into 



a continuum of surface-projected bulk band states in the limit 

d — !■ oo. 

The question arises, when do we have to treat the electron 
states as confined in a slab and when as bulk states. The deci- 
sive quantity is the electron coherence length £c- If (-c "C d, it is 
equivalent to the bulk mean free path and is determined by bulk 
scattering centres. When the coherence length exceeds the slab 
thickness, Ic > d, the electrons between the two slab bound- 
aries behave as standing waves in the normal direction and as 
extended states in the other two directions. Therefore, in case 
of ultrathin slabs, one has to distinguish between the electron 
coherence length in the direction normal to the film, and the 
electron mean free path in the plane of the film, imfp, which is 
predominantly limited by the surface roughness. 

The DFT calculations were performed on 1 to 8 monolay- 
ers thick, (lll)-oriented Pb slabs, separated by 10 A of 
vacuum, and with periodic boundary conditions. If the vac- 
uum layer were too thin, we would observe dispersion of the 
electron bands in the z direction. Indeed, we observed no dis- 
persion in the occupied bands. The electron band energies and 
the total energy were calculated ab-initio with the full-potential 
linearised augmented plane-wave method in the local-density 
approximation [19] as implemented in the WIEN2k code [20]. 
A mixed basis set of augmented plane waves plus local orbitals 
(APWh-1o) [21] for low orbital momenta (/ < 2) and linearised 
augmented plane waves (LAPW) for all the higher orbital mo- 
menta were used. The spin-orbit interaction was included self- 
consistently by applying the second- variational method and us- 
ing the scalar-relativistic eigenfunctions as basis [22]. The Pb 
muffin-tin radius was set to 2.6 a.u. and a tetrahedral mesh 
of about 600 fc— points in the irreducible part of the Brillouin 
zone. Fig. 1, was used in the self-consistent electronic structure 
calculations. The kinetic-energy cutoff was set to ii^,^ax = 9.5 
Ry and the plane-wave expansion cutoff to = 196 Ry. 

Later, in the calculation of Fermi lines and conductivities, we 
used a mesh with about 4000 k points in the irreducible part of 
the Brillouin zone. 

First, the energy of bulk Pb FCC crystal was minimised 
to find the equilibrium bulk lattice constant, Ocaic = 4.88 A, 
which is 1 .4% smaller than the experimental room-temperature 
value Ooxp = 4.95 A and 3% smaller than the calculated value 
reported by Materzanini et al. who used the Cambridge Serial 
Total Energy Package [15]. For the slab calculations, we used 
the corresponding in-plane hexagonal lattice constant, a = acaic/V2 ^ 
3.45 A. Except for d = 1 ML, no visible change in the band- 
structure is observed if either the experimental lattice constant 
was used in the simulations or if the lattice was unrelaxed. 

The transport properties are calculated using the Boltzmann 
transport equation for two-dimensional systems in the relaxation- 
time approximation. The two-dimensional description is nec- 
essary when the charge-carrier mean free path exceeds the slab 
thickness. As we shall see later, this condition is fulfilled for 
all but the thinnest (d = 1) slabs. For polycrystalUne or amor- 
phous metallic films the transport properties are isotropic in the 
plane of the surface and the conductivity is [23] 
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where e is the electron charge, v„ the group velocity of an elec- 
tron in the subband n at the Fermi energy ep and T„(eF) its 
relaxation time. The sum runs over all subbands crossing the 
Fermi energy and the integral is along each Fermi line. The 
spin degeneracy is included in the prefactors. The Hall conduc- 
tivity is [23] 

The mechanisms influencing the relaxation time are elastic 
scattering on bulk imperfections and elastic scattering on rough 
surfaces, l/r = l/rf, -I- 1/ts- Inelastic scattering on phonons 
can be neglected if the conductivities are investigated at low 
temperatures. In ultrathin films, the main scattering mechanism 
is elastic scattering on surface roughness for which r„(eir) is 
given by 

= ¥ E \{n'k'\ni\nk)fS{en'{k') - ep) (3) 

where {n'k'\T-Li\nk) is the scattering matrix element. For small, 
uncorrelated islands, the matrix element is independent of n 
and k, it is equal to U'^pg/ Sd"^ where U is the strength and pg 
the surface density of the scattering centres [10]. S* is the area 
of the film. If this situation is fulfilled, the scattering is diffuse 
[4] and the scattering rate simplifies to 

1 ^ u^pg ^ r 

T 2nh?d^ ^ f \vnieF)\- 

In this case the charge-carrier momentum is lost completely, 
scattering causes transitions between all subbands that cross 
€f and also the relaxation time of the electrons at is inde- 
pendent of the subband index n and electron momentum k. In 
general, vp is only weakly /c— dependent so that r ex d^/L re- 
sults, where L is the length of all Fermi Unes. As we shall see 
later, L and therefore also r and cto are approximately Unear 
functions of d. 



3 Results 

We start our analysis by presenting the results of the lattice 
relaxation and bandstructure calculations of relaxed Pb(lll) 
slabs. The thickness dependence of the interlayer spacings is 
shown in Fig. 2. The topmost layer relaxes always inwards. In 
addition, we find a quantum-size induced oscillation of dij as 
a function of thickness; we observe simultaneous inward re- 
laxation of all investigated layers for d = 2,5 and 7 ML thick 
films. Whereas the relaxations of the topmost layers agree quite 
well with the values reported in Ref. [15], the relaxations of 
the inner layers differ. Our spacings ^23 show similar QSE os- 
cillations as in Ref. [15], but they are systematically smaller 
(~ 2% of dbuik, c?buik is the bulk interlayer spacing in the (111) 
direction). The amplitude of the QSE oscillations in ^34 is ^ 
1.5% of dbuik and is thus comparable to the amphtude of ds4 



in Ref. [15], but with opposite phases as a function of layer 
thickness. 

Fig. 3 shows the excess slab energies Sex = £^siab— c^-^buik, 
where EsUh is the total energy of a slab and i?buik = —41834.809 
Ry the energy/atom in the bulk. Clearly seen are the QSE- 
induced oscillations in with minima for 2, 4 and presum- 
ably also 8 ML thick slabs, in partial agreement with the DFT 
calculations of Ref. [15]. 

The electron energy bands along the main symmetry di- 
rections of the two-dimensional hexagonal BriUouin zone for 
d = 4 are shown in Fig. 4. Each band is split into d subbands 
which do not cross because of spin-orbit interaction. The num- 
ber of subbands crossing e^? is roughly proportional to d. This 
property will have direct consequences on r and on the con- 
ductivities. The points of subband crossings with cp form two- 
dimensional closed loops, i.e., Fermi lines. Depending on the 
occupancy of the states enclosed by a Fermi line we have hole- 
or electron-like carriers. Their ratio will play a crucial role in 
determining the sign of the Hall coefficient Rh = (th/ctq. 
Also shown in Fig. 4 is the electron density of states. Close to 
the bottom of the subbands the density of states has steps, spe- 
cific for parabolic bands in two dimensions. Proportional to d is 
also the total length of Fermi lines, L, in the two-dimensional 
BriUouin zone, see Fig. 5. This, together with diffuse scatter- 
ing, is the main reason why ao is approximately proportional 
to d. 

We next calculated the relaxation time, Eq. 4, assuming a 
potential depth 3 eV, a terrace height one ML (- 3 A), a terrace 
area 100 (U = 900 eVA^) and a terrace density pg = 10""' 
A~^, see Fig. 6. These numbers were chosen in order to sim- 
ulate closely the experimental situation described in Ref. [24] 
and to obtain a close fit to the experimental values for ctq and 
Rh [24]. We note that elastic scattering at the surfaces causes 
transitions between all subbands that cross the Fermi energy. 
The summation over all final states in the scattering rate (4) 
causes t cx d. As we see from Fig. 6, r cx d — 1 describes the 
thickness dependence better. The reason is found in the details 
of the bandstructure (see below). The deviations from Unearity 
at d = 3, 6 and 8 ML appear when the Fermi energy coincides 
with a local peak in the electron density of states, thus opening 
more channels for elastic scattering. Also shown in Fig. 6 is the 
"experimental" value t'^'P, 

^exp ^ 47r^;id gg^P 

§ dl\vn{eF)\ ' 

obtained from Eq. (1) by combining a'^'^ from Ref. [18] with 
the calculated value for ^ dl\vn{eF)\- We also calculated ao, 
Fig. 7, and (Th, Fig. 8. (Tq has a very similar thickness depen- 
dence as T, indicating, that r is the main source of the thickness 
dependence of ctq. 

The Hall conductivity, an, on the other hand, is a compU- 
cated function of d, mainly because of compensating effect of 
electron and hole subbands contributions. 

There have been several attempts to relate the Hall coeffi- 
cient to the electron density of states, in particular for disor- 
dered metals [25]. Therefore we decided to test the above rela- 
tions. In Fig. 9 we show the dependence of calculated tJo/r 
and (Tjj/t^ on the position of the Fermi level for a d = 4 
ML slab. We artificially shift the Fermi energy and calculate 
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the corresponding ai^/T{eF) and (Jh /T^i^p)- As is seen from 
this figure, ctq /r is a relatively smooth function of whereas 
ohIt^ oscillates strongly with ep and changes sign: it is neg- 
ative at the bottom of each subband, when the charge carri- 
ers are electron-like, and positive at the top of the subbands 
when the carriers are hole-like. In between, they are compli- 
cated functions of the Fermi energy and reflect different singu- 
larities in the bandstructure. As a consequence, the dependence 
of Rh on is given predominantly by gh /t^ ■ Comparison of 
fff/T^ (ef ). Fig. 9, and of the density of states. Fig. 4(b), shows 
that for d = 4 ML, there is no evident relation between Rh and 
the density of states. Such a relation would be valid only if we 
had a single three-dimensional band, i.e. when the charge car- 
riers at the bottom (top) of the band are electrons with Rh < 
(holes with Rh > 0) and at the same time the density of states 
is a rising (falhng) function of the energy. In case of several 
inter-penetrating bands, as is the case of thin films, such a sim- 
ple picture breaks down [26-28], as is also obvious from our 
results. We also find no obvious relation between cfq and the 
density of states. 

4 Comparison with the experiments 

In Figs. 6 and 7 we compare calculated t and (Tq with the exper- 
iments. It has been observed experimentally that the transport 
properties of Pb films are almost independent of the substrate 
reconstruction [18]. Therefore, metals on Si can be considered 
to a reasonable first approximation as free-standing metal films. 
This is not true for the first Pb layer, which will be discussed 
separately below. 

Clearly seen in theory and experiments is an approximately 
linear thickness dependence of (To which corroborates the pic- 
ture that the scattering is dominated by surface roughness and 
that the roughness does not change with film thickness. 

We identify several conductivity regions. Below the per- 
colation threshold at d ~ 0.8 ML the charge carriers are lo- 
calised, the transport is strongly temperature activated and ctq 
is so small that it is of no relevance for our considerations. 

Between the percolation threshold and 1 ML the Pb films 
are pseudomorphous, i.e. the lattice constant is ^ 10% larger 
than the corresponding bulk Pb value. We performed also DFT 
and conductivity calculations for a d = 1 ML thick Pb slab 
with the lateral lattice constant of a Si(lll) substrate and we 
find that the electrical properties of a 1 ML film are extremely 
sensitive to the lateral lattice constant a. A 10% increase in 
the lattice constant changes the metal to an indirect bandgap 
semiconductor. The experimental conductivity of a 1 ML film 
is very small and has an approximately linear temperature de- 
pendence [18], supporting the theoretical conjecture that the 1 
ML Pb(lll) films are semiconducting. 

Between ~ 1 ML and 3.5 ML, the annealed films are 
amorphous and CTq'^p exhibits a small, almost linear increase 
with temperature up to ^ 100 K, indicating a small activated 
contribution. Films with d > 4 ML are crystalUne and a'^'^ de- 
creases with increasing temperature (phonons). We could not 
observe any influence of disorder on the conductivity of ultra- 
thin films. The experimentally observed dip in cto at 3 ML (see 
Fig. 7) is fully reproduced in our calculation with ideally or- 
dered slabs. It is related to the electronic bandstructure, but has 



nothing to do with disorder in the film. Therefore, it is a pure 
quantum-size effect. 

The thickness dependence of ctq is, like r, oc d — 1. It is 
governed by scattering on the surface or interface roughness - 
bulk disorder would give a thickness-independent ctq . The main 
source of the oscillations of ctq with d is the quantisation of the 
Pb band structure in the direction normal to the film together 
with the strong spin-orbit interaction. 

Our results for o-q are consistent with the picture that the 
electrons are subject to diffuse elastic scattering at the sur- 
faces of the films [1,2,4], in agreement with the original Fuchs- 
Sondheimer model, but they are at variance with the model 
of scattering at surfaces with uncorrected roughness, where 
a oc (P was predicted [10]. 

The calculated Hall conductivity is compared to cr^'' in 
Fig. 8. The agreement is good except for d = 5 ML where the 
discrepancy most probably lies in the lifetime broadening of 
the electron states [18] which is neglected in the present analy- 
sis. 



5 Discussion 

Our calculation used for the interpretation of the experimental 
data is based on a DFT band structure calculation combined 
with the Boltzmann transport equation. The only adjustable 
parameter, used in the calculation, is the scattering strength 
U'^Ps- The physical scenario is described by electronic states 
in ultrathin films that are localised in the direction normal to 
the surface and are scattered by surface roughness. The elec- 
tron states form discrete subbands since the electron coherence 
length, limited by scattering in the film bulk, exceeds the film 
thickness. 

The results of the calculation are consistent with the pic- 
ture that the electrons are subject to diffuse elastic scattering 
on film surfaces [1, 2,4] for which ctq and r are proportional to 
the number of conducting metallic layers, in our case {d— 1), in 
agreement with the original Fuchs-Sondheimer model and with 
the experiments [18]. The (lateral) mean free path of electrons 
in Pb ultrathin films on Si(l 1 1) exceeds the Pb bulk lattice con- 
stant when d > 1 ML. According to the loffe-Regel criterion 
[29], ^ ~ 1/A:f ~ flfc (fflb is the bulk lattice constant), the elec- 
trons in d > 1 ML thick films are not localised, and the use of 
the Boltzmann transport equation is justified. 

The essential difference between the Fuchs-Sondheimer and 
the present model is that in the former the electrons are de- 
scribed by three-dimensional Bloch states whereas in our ap- 
proach the electron states are described by two-dimensional 
momenta and discrete subbands. The former approach is 
adequate for thick slabs, whereas our approach is valid for ul- 
trathin slabs when the electron coherence length exceeds the 
slab thickness d. 

In experiment, no free-standing films can be investigated, 
i.e. an insulating substrate is always required that tries to im- 
pose its own periodicity on the adsorbed film, depending on the 
interface energies, and on lattice mismatch. For a perfect film, 
an atomically sharp interface and a well defined periodicity at 
the interface are required (e.g., in form of a superstructure). 
In the Pb/Si(ll 1) case, where the lattice mismatch is approxi- 
mately 10%, no lattice match, apart from the first monolayer. 
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can be expected. This is the reason for the formation of the 
amorphous films up to 4 monolayers, whose many degrees of 
freedom are obviously needed to allow for the growth of a Pb 
film with its own lattice constant for thicker layers. Neverthe- 
less, it seems that the short range order in these amorphous 
films is still sufficiently good to allow a comparison of conduc- 
tivity and Hall conductivity calculated for perfect crystalline 
order with experimental data obtained in amorphous films. The 
calculations with perfectly ordered films even reproduce details 
of the conductivity as a function of layer thickness. 

As seen from our analysis, the diffuse interface scatter- 
ing mechanism still limits the conductivity, even for the crys- 
talline films at larger d. For crystalline layers the interface to 
the vacuum is atomically smooth for complete layers, but can 
be roughened by adsorption of incomplete layers, as seen by a 
reduction in conductivity [24] . These roughness-induced changes 
of conductivity, however, are small. Therefore, the Pb-vacuum 
interface does not seem to provide the main contribution to dif- 
fuse scattering for the crystalline Pb films (the terraces are too 
large and their density too small). An electronically still very 
rough interface, on the other hand, may exist between the crys- 
talline and the amorphous layers. Because of the lack of or- 
der in the amorphous state, this interface scatters at random in 
all directions, i.e. it still acts as diffusely scattering interface. 
Additional contributions to electron scattering from the Si-Pb 
interface are conceivable. 

The single pseudomorphous monolayer of Pb on Si(lll) 
has a narrow indirect band gap and is thus semiconducting, ac- 
cording to the DFT calculations, in agreement with the experi- 
ments, where the conductivity has an approximately linear tem- 
perature dependence [18]. This is not the case for free-standing 
Pb monolayers, which correspond to Pb deposited on semicon- 
ducting 1 ML Ag-covered Si(lll)(\/3 x ^/3) [17], where Pb 
still has a Fermi Une and is thus, strictly speaking, a metal, but 
with a very small conductivity. 

6 Conclusions 

Our analysis of the charge-carrier scattering and transport in 
ultrathin metallic Pb films has shown that they are very sensi- 
tive to the reduced dimensionality of ultrathin layers resulting 
in classical and quantum size effects. The relaxation time r 
and the conductivity ctq are Umited by the scattering on rough 
surfaces and are approximately linear functions of d (classi- 
cal size effect) with quantum-size induced deviations from lin- 
earity at d = 3, 6 and 8 ML. The proportionaUty constant in 
fjo oc {d — 1) depends on roughness of the interface to the 
substrate. The deviations for d = 3,6 and 8 ML appear when 
the Fermi energy coincides with a local maximum in the elec- 
tron density of states and causes a drop in r. These deviations 
are the analogue of the condition id « j^^Ff^ (where i and j 
are integers and the Fermi wave length) in the free electron 
model [5,6]. Rh is very sensitive to details of the bandstruc- 
ture and to a large extent insensitive to disorder in the film and 
to surface roughness. 

As demonstrated in this paper, the charge-carrier transport 
is very suitable for investigating the QSE, since it is to a large 
extent insensitive to disorder in the film and is limited by the 
interface roughness. 



The quantum-size-induced oscillations of aH and Rh are 
strong and originate from partially compensating contributions 
to the Hall conductivity that emerge from different subbands 
crossing the Fermi energy. The advantage of Rh over cth is 
that Rh is almost insensitive to the charge-carrier lifetime and 
consequently to ordering in the film. 

Whereas one ML of free-standing Pb slab is metallic, a sin- 
gle monolayer of Pb, pseudomorphous with Si(lll), is semi- 
conducting. In both cases, however, the conductivity and the 
mean free path are small. The band gap in the pseudomorphous 
case is caused by spin-orbit splitting. 

The method applied here can be used for other ultrathin 
metallic films, provided the electron coherence length (or mean 
free path) is larger than the film thickness. We have shown 
that the quantum-size effects depend on the details of the elec- 
tronic bandstructure and cannot be treated quantitatively with 
the free electron models, in particular in case of heavier metals 
when the spin-orbit interaction is strong. This applies also to 
aH which is very sensitive to details of the bandstructure. 

This work was supported by grants from the Deutsche Forschungs- 
gemeinschaft, from the Deutsche Akademische Austauschdienst and 
from the Slovenian MSZS. 
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Fig. 3. Thickness dependence of the excess slab energies Eex of re- 
laxed symmetric Pb(l 1 1) slabs. 
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Fig. 1. Irreducible part (dark grey) of the two-dimensional hexagonal 
first Brillouin zone of a free-standing slab. For the slab on a substrate 
the in-plane mirror symmetry is broken and the irreducible Brillouin 
zone comprises both grey areas. 



Fig. 4. Calculated bandstructure (left) and density of states (right) of a 
free-standing 4 ML thick Pb(l 11) film. The shaded areas in the band- 
structure represent the surface-projected bulk electron energy bands, 
whereas the shaded areas in the density-of-states are the occupied 

states of the 4 ML thick Pb film. 




Fig. 2. Interlayer spacings dij (in A) of relaxed symmetric Pb(lll) 
slabs as a function of slab thickness. The lowest indices correspond to 
the surface layers. Dashed line indicates the relaxed bulk spacing. 




Fig. 5. Thickness dependence of the total length L of calculated Fermi 
Unes in the two-dimensional Brillouin zone. Dashed Une is a Unear 
fit through the the calculated points. Also shown is the calculated 
Fermi-line length of a d = 1 ML thick film, pseudomorphous with 
theSi(lll) substrate. 
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Fig. 6. Thickness dependence of the charge-carrier scattering relax- 
ation time r assuming equal roughness (equal U^ps) of all slab thick- 
nesses d. In case of diffuse elastic scattering, t is independent of the 
subband index n and electron momentiun k. Experimental points are 
given by Eq. (5) where the experimental conductivities are from 
Ref. [18]. Dashed line is a linear fit to the experimental points. 




d (ML) 



Fig. 7. Thickness dependence of the calculated [Eq. (1)] and low- 
temperature experimental [18] conductivity ctq. Notice the remarkable 
similarity between ao and r which is the consequence of linear thick- 
ness dependence of L, Fig. 6. Dashed line is a linear fit through the 
experimental points. (In Fig. 2 of Ref. [18] ao must be multiplied by 
a factor 100.) 
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Fig. 8. The calculated Hall conductivity an, Eq. (2), is compared to 
the experimental Hall conductivity, ctJJ'' = Ruicrf^^)^ . Huge oscilla- 
tions are the effect of cancellations between electron and hole contri- 
butions. 
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Fig. 9. Energy dependence of the calculated electrical conductivity 
ao/r and Hall conductivity cth / for a 4 ML thick slab. 



